Quantum optics with quantum gases: 
controlled state reduction by designed light scattering 
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Cavity enhanced light scattering off an ultracold gas in an optical lattice constitutes a quantum 
measurement with a controllable form of the measurement back-action. Time-resolved counting 
of scattered photons alters the state of the atoms without particle loss implementing a quantum 
nondemolition (QND) measurement. The conditional dynamics is given by the interplay between 
photodetection events (quantum jumps) and no-count processes. The class of emerging atomic 
many-body states can be chosen via the optical geometry and light frequencies. Light detection 
along the angle of a diffraction maximum (Bragg angle) creates an atom-number squeezed state, 
while light detection at diffraction minima leads to the macroscopic superposition states (Schrodinger 
cat states) of different atom numbers in the cavity mode. A measurement of the cavity transmission 
intensity can lead to atom-number squeezed or macroscopic superposition states depending on its 
outcome. We analyze the robustness of the superposition with respect to missed counts and find 
that a transmission measurement yields more robust and controllable superposition states than the 
ones obtained by scattering at a diffraction minimum. 

PACS numbers: 03.75.Lm, 42.50.-p, 05.30.Jp 
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I. INTRODUCTION 

Both, quantum optics and physics of ultracold quan- 
tum gases represent nowadays the well-established and 
actively developing fields of modern quantum science 
d 0. However, the interaction between the two fields 
is far from being close. 

Historically, classical optics treating the light as clas- 
sical electromagnetic waves has become one of the most 
developed and fruitful fields of physics. It has provided 
us a lot of technological breakthroughs, e.g., the highest 
level of measurement precision. Quantum optics, which 
considers the light as quantum particles (photons), thus 
going beyond the mean-field (classical) description of the 
electromagnetic waves, currently is also a well-developed 
field, both theoretically and experimentally [![. 

The progress in laser cooling techniques in the last 
decades of the 20th century led to the foundation of a new 
field of atom physics: atom optics. It was shown that the 
matter waves of ultracold atoms can be treated similar 
to light waves in classical optics and can be manipulated 
using the forces and potentials of laser light beams. The 
quantum properties of matter waves beyond the mean- 
field description became accessible after 1995, when the 
first Bose-Einstein condensate (BEC) and many other 
fascination quantum states of bosonic and fermionic ul- 
tracold atoms were obtained [1 Q ■ An exciting demon- 
stration of "quantum atom optics" was presented in 2002, 
when the quantum phase transition between two states of 
atoms with nearly the same mean density, but radically 
different quantum fluctuations was obtained:_superfluid 
(SF) to Mott insulator (MI) state transition 
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The roles of light and matter in optics and atom optics 
are completely reversed. Various devices known in optics 
as beam-splitters, mirrors, diffraction grating, etc. are 
created using light forces and applied for matter waves. 
However, up to now, the absolute majority of even very 
involved setups and theoretical models of quantum atom 
optics treat light as an essentially classical axillary tool 
to prepare and probe intriguing atomic states. In this 
contexts, the periodic micropotentials of light (optical 
lattices) play the role of cavities in optics enabling one 
to store and manipulate various atomic quantum states. 

Quantum optics with quantum gases should close the 
gap between quantum optics and atom optics by address- 
ing phenomena, where the quantum natures of both light 
and matter play equally important role. Experimentally, 
such an ultimate quantum level of the light-matter in- 
teraction became feasible only recently, when the quan- 
tum gas was coupled to the mode of a high-Q cavity 
d, 0, H, Q ■ Even early theoretical works on scattering 
of quan tized li ght from a BEC was not realized so far 
0, lH E EI 13, EE M El El, El, M SO, HI . How- 
ever, it is cavity quantum electrodynamics (QED) with 
quantum gases that will provide the best interplay be- 
tween the atom- and light-stimulated quantum effects. 

On the one hand, the quantum properties of atoms 
will manifest themselves in the scattered light, which 
will lead to novel nondestructive methods of probing 
and manipulating atomic states by light measurement 

di m, si si m m, m m m m m, bhi. o n 

the other hand, the quantization of light (i.e. trapping 
potentials) will modify atomic manybody dynamics well- 
known only for classi cal p otentials and give rise to novel 
quantum phases [H, M, H H, EE El El El ■ This 
paper addresses the first problem. 

In this work, we will extend our treatment of the de- 
tection of light scattered from ultracold atoms in opti- 



cal lattices presented in Ref. |26|], where the quantum 
measurement of light was considered as a method to pre- 
pare particular quantum states of atoms thanks to the 
measurement back-action. Except detailing the theoret- 
ical approach, we will consider a different optical con- 
figuration of the measurement process: cavity transmis- 
sion spectroscopy instead of light scattering. As we will 
show, it turns out, that such a configuration will allow us 
much better control and flexibility in atomic state prepa- 
ration. This is in particular true for the preparation of 
macroscopic superposition states (known as Schrodingcr 
cat states), which can be more robust in comparison to 
those considered before [H, E3 • 

We will show that detecting the light scattered from ul- 
tracold atoms in an optical lattice enables one to prepare 
various types of the atom-number squeezed and macro- 
scopic superposition states. An important point, is that 
the type of the state and its properties strongly depend 
on the optical geometry. Varying the optical parame- 
ters (angles between the laser beams and lattice atoms, 
light frequencies, or lattice period) one can prepare vari- 
ous quantum states of ultracold atoms. Moreover, as the 
optical measurement is nondestructive, in the sense of 
the quantum nondemolition measurements (QND), one 
can make sequential measurements on the same sample 
without completely destroying its quantum properties. 

Generalizing the methods developed for spin squeez- 
ing in thermal atomic ensembles Q EE EE E3, EH for 
the case of light scattering from ultracold quantum gases 
will enable ones to reach the unprecedented level of the 
measurement precision, which is required, e.g., for atomic 
clocks and the detection of gravitational waves. 

The paper is organized as follows. In Sec. II, a gen- 
eral theoretical model is formulated, which, in Sec. Ill, 
is extended to the quantum back-action of the light mea- 
surement on the atomic state. Section IV is devoted to 
the scattering of the transverse probe in to the cavity. In 
Sec. V, the properties of photon statistics are analyzed. 
The state preparation by the cavity transmission mea- 
surement is discussed in Sec. VI. The robustness of the 
macroscopic superposition states is analyzed in Sec. VII. 
The main results are summarized in Sec. VIII. 



II. GENERAL MODEL 

We consider N ultracold atoms trapped in a periodic 
optical lattice potential with M sites. In addition to the 
strong trapping beams, a subset of K sites, with K < M, 
is illuminated by a weak probe at some angle to the lat- 
tice atoms. The scattered light is collected by a cavity. 
The photons leaking out of the cavity are counted then 
with a photodetector. The schematic setup is shown in 
Fig. 1, where, for simplicity, the trapping beams are not 
shown. In general, the setup is not restricted by a sim- 
plified scheme in Fig. 1: the lattice can be one-, two-, 
or three-dimensional; the probe and cavity modes can be 
standing or traveling waves, or they can be even formed 
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FIG. 1: (Color online) Setup. A lattice is illuminated by 
the transverse probe ao and probe through a mirror 77. The 
photodetector measures photons leaking the cavity. Due to 
the quantum back-action, the light measurement leads to the 
modification of the atomic quantum state. 



by two different modes of the same cavity. K illuminated 
sites can be selected in some nontrivial way, e.g., each sec- 
ond lattice site can be illuminated by choosing the probe 
wavelength two times larger than the lattice period. In 
addition to the external probe, a probing through a cav- 
ity mirror will be considered. 

The main goal of this paper is to study how the mea- 
surement of the photons leaking out of the cavity will 
affect the quantum state of the ultracold atoms. Such 
a measurement back-action is possible due to the entan- 
glement between light and matter, which develops during 
the interaction process. According to quantum mechan- 
ics, in the presence of the entanglement, the measurement 
of one of the quantum subsystems (light) will also affect 
another quantum subsystem (atoms). 

We will use the physical model presented in details in 
Ref. . The theoretical formulation starts with a gen- 
eralized Bose-Hubbard model including the quantization 
of light. As here, unlike Ref. [4!|, the far off-resonant 
light scattering will be considered, the role of the atomic 
excited states is not important and they can be adia- 
batically eliminated from the dynamics. As we will be 
focused here on the quantum measurement process, it 
is reasonable to neglect several processes have been con- 
sidered in other works [33, E3, EH to avoid the complica- 
tions. One of such processes to be neglected is the details 
of atomic dynamics. The tunneling of atoms in the lattice 
potential plays indeed an important role in establishing 
a particular quantum atomic state. However, after the 
state is established, one can assume that the scattering 
of the probe occurs on the time-scale faster than slow 
tunneling. Thus, from the light scattering point of view, 
the atomic distribution can be considered frozen and the 
tunneling is not important. It is especially reasonable to 
neglect the tunneling dynamics in this paper, because we 
will show that even after that rather obvious dynamics 
is neglected, there is still non-trivial dynamics (quantum 
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jumps) exclusively associated with the quantum measure- 
ment process, which is the main subject of our present 
work. 

After those simplifications, the Hamiltonian of the 
problem takes the form: 

H = %{u)x + UnDu)a\ai + ?L£/io(-D 10 ao a i + D w a a\) 

-ih(rj*ai - r?aj),(l) 

where a\ is the cavity-mode annihilation operator of the 
frequency u>\ . The external probe of the frequency lo v is 
assumed to be in a coherent state, thus its amplitude is 
given by a c- number ao- The spatial mode functions of 
the probe and cavity modes are ui,o(r). Ui m = gig m /A a 
(I, tyl — 0, 1), where 51,0 are the atom-light coupling con- 
stants, A a = lu\ — uj a is the cavity-atom detuning, rj 
is the amplitude of the additional probing through a 

' p . We assumed the probe- 
— u>\ <C A a . The operators 
u sum contributions from all 



mirror at the frequency w 
cavity detuning A p = u> 

Dim = EjLl_<( r iK . 

illuminated sites with the atom-number operators hj at 
the position rj . 

The first term in Eq. |T]) describes the atom-induced 
shift of the cavity resonance. The second one reflects 
scattering (diffraction) of the probe ao into a cavity mode 
ax- For a quantum gas the frequency shift and probe- 
cavity coupling coefficient are operators, which leads to 
different light scattering from various atomic quantum 
states [H 0, 12II23. 

The Hamiltonian {T]) describes QND measurements of 
the variables related to Di m measuring the photon num- 
ber ajai, as the criteria for the QND measurements are 
fulfilled [50]. Note, that one has a QND access to vari- 
ous many-body variables, as Di m strongly depend on the 
lattice and light geometry via uo,i(r). This is an advan- 
tage of the lattice comparing to single- or double-well 
setups, where the photon measurement back-action was 
considered [UH^H^Hil]. Moreover, such a geometrical 
approach can be extended to other quantum arrays, e.g. 
ion strings 
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For example, one can consider a ID lattice of the pe- 
riod d with atoms trapped at Xj = jd (j = 1, 2, .., M). In 
this case, the geometric mode functions can be expressed 
as follows: Uo t i(rj) = enp(ijko t i x d + <j>o,ij) f° r traveling 
waves, and ito,i(rj) = cos(jko.i x d + (/>o,ij) f° r standing 
waves, where feo,ix = |ko,i | sin #0,1, #0,1 are the angles 
between mode wave vectors ko,i and a vector normal to 
the lattice axis. In the plane-wave approximation, addi- 
tional phases 4>o,ij are j-independent. The general angu- 
lar distributions of light scattered from ultracold atoms 
in optical lattices were presented in Refs. [13, HH, H3] 

For some geometries, Dn can reduce to the operator 
Nk = J2f=i ™j 01 tne atom number at K sites [23I.I24I.I25I] 
(if a± is a traveling wave at an arbitrary angle to the lat- 
tice, or the standing wave with atoms trapped at the 
antinodes). If the probe and cavity modes are coupled 
at a diffraction maximum (Bragg angle), i.e., all atoms 



scatter light in phase, u*(rj)uo(rj) = 1, the probe-cavity 

coupling is maximized, D\q = Nk- If they are coupled at 
a diffraction minimum, i.e., the neighboring atoms scat- 
ter out of phase, L>io = Sj=i ( — is the operator 
of number difference between odd and even sites. Thus, 
the atom number as well as number difference can be 
nondestructively measured. Note, that those are just two 
of many examples of how a QND-variable, and thus the 
projected state, can be chosen by the geometry. 



III. MEASUREMENT BACK-ACTION 

In this section we present the solution for the quantum 
state of the coupled light-matter system including the 
measurement process. Importantly, it is possible to ob- 
tain an analytical solution with a very transparent phys- 
ical meaning thanks to the approximations used (slow 
tunneling and coherent state of the external probes). 

The initial motional state of the ultracold atoms 
trapped in the periodic lattice potential at the time mo- 
ment t = can be represented as 



l* a (0)> 



E 



c° q \qu 



,QM) 



(2) 



which is a quantum superposition of the Fock states cor- 
responding to all possible classical configurations q = 
{qi, .., 5m} of N atoms at M sites, where qj is the atom 
number at the site j. For each classical configuration q, 
the total atom number is conserved: qj = N. This 
superposition displays the uncertainty principle, stating 
that at ultralow temperatures even a single atom can be 
delocalized in space, i.e., there is a probability to find an 
atom at any lattice site. We will show, how this atomic 
uncertainty is influenced by the light detection. 

For example, for a limiting case of the MI state, 
where the atom numbers at each lattice site are pre- 
cisely known, only one Fock state will exist in Eq. ©: 
I^mi) = 1 1 , 1 , - - , 1) for the MI with one atom at each 
site. On the other hand, the SF state is given by the 
superposition of all possible classical configurations with 
multinomial coefficients: 



l*SF> = 



(VM) 



N 



E 



AT! 



qik2 l - 



\qi,q 2 ,--qM) 



(3) 



Thus, the atom number at a single site as well as the 
atom number at K < M sites are uncertain in the SF 
state. 

As an initial condition, we assume that at the time 
moment t = the light and matter are disentangled, and 
the initial state of light is a coherent state with the am- 
plitude ao- Thus, the initial quantum state of the system 
is given by the product state |*(0)) = |* a (0))|a ). In 
particular, initially, the light can be in the vacuum state 
|0). 
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We use the open system approach [56| to describe the 
continuous counting of the photons leaking out the cavity 
of the cavity decay rate k. According to that approach, 
when the photon is detected at the moment U, the quan- 
tum jump occurs, and the state instantaneously changes 
to a new one obtained by applying the cavity photon an- 
nihilation operator |\& c (tj)) — > ai\^B c (ti)) and renormal- 
ization (the subscript c underlines that we deal with the 
state conditioned on the photocount event). Between the 
photocounts, the system evolves with a non-Hermitian 
Hamiltonian H — i%Ka\a\. Such an evolution gives a 
quantum trajectory for \^ c (t)) conditioned on the detec- 
tion of photons at times ti,t 2 , .... The probability of the 
photon escape within the time interval t is 2Kt{a\a\) c , 
where {a\a\) c is the conditional photon number in the 
cavity, i.e., the photon number calculated for the condi- 
tional quantum state Y$> c (t)). 

The state \^ c {t)) should be found by solving the 
Schrodinger equation with the non-Hermitian Hamilto- 
nian for no-count intervals and applying the jump op- 
erator a\ at the moments of photocounts. Thanks to 
the slow tunneling approximation, the Hamiltonian ([T]) 
does not mix the Fock states in the expression (|2|). So, 
the problem is significantly reduced to separate find- 
ing solutions for each classical atomic configuration q = 
{qi, ..,qu}, after that the full solution will be given by 
the superposition of those solutions. 

The next simplification appears thanks to the use 
of the external probes in the coherent state. It is 
known [5^ . l58j that, if a coherent probe illuminates a 
classical atomic configuration q in a cavity, the light re- 
mains in a coherent state, however, with some pre-factor: 
exp[<i> g (?j)]|a ? (?j)). The pre-factor is indeed not impor- 
tant for a single classical configuration as it disappears 
after the renormalization, but it will play a role, when 
the superposition of classical solutions with different pre- 
factors will be considered. Moreover, the light ampli- 
tude a q (t) is simply given by the solution of a classical 
Maxwell's equation: 



a q {t) 



fj - iUi a D q 10 
i(U n D g n - A p ) + 



_ e -iuJ P t _|_ 



«0 



fj - iU w a a D\ { 
i{U u D g n - A p ) 



(4) 



where we introduced the constant probe amplitudes ao 
and fj as ao = fioexp(— iuj p t) and rj = ijexp(—iu p t); 
Dim — 2j=i u i{ v j) u m{^j)qj is a realization of the oper- 
ator Di m at the configuration q = {gi, ..qM }■ As in classi- 
cal optics, the first term in Eq. ([4]) gives the oscillations at 
the probe frequency, while the second term gives the tran- 
sient process with the oscillations at the cavity frequency 
shifted by the dispersion, ui + UnD^, which decays with 
the rate k. In the following, we introduce the slowly vary- 
ing light amplitude a q (t) as a q (t) — a q (t) exp (—iuj p t) 
and, for the notation simplicity, will drop the tilde sign 
in all amplitudes ao, fj, and a q (t). 



The function $ q (i) in the pre-factor is a complex one 
and is given by 



l(w* q 



iU 10 a D q 1Q a* o 



c.c. 



n\a q \ 2 



dt,(5) 



where the light amplitude a q (t) is given by Eq. |4|. 

First, let us consider the solution for the atomic state 
initially containing a single Fock state |gi, .., gjvf)- The 
solution for light is given by the solution for the classi- 
cal configuration q = {gi, .., 5m}- So, the evolution is 
given by the product state |gi, .., qM)\oi q (t)). An impor- 
tant property of this solution is that a quantum jump 
does not change the state, since applying the jump oper- 
ator ai simply leads to the pre-factor a q (t), which disap- 
pears after the renormalization. Therefore, even in the 
presence of photocounts (i.e. quantum jumps), the time 
evolution is continuous and is given by Eq. ((4]): after a 
transient process for t < 1/k, the steady state for a q (t) 
is achieved. Note, that in contrast to many problems 
in quantum optics [561 ]. where the steady state is a re- 
sult of averaging over many quantum trajectory, here, 
the steady state appears even at a single quantum tra- 
jectory for t > 1/k. This is a particular property of 
the coherent quantum state. The continuity of evolu- 
tion of the state |gi, .., qM)\ot q (t)), and, hence, more gen- 
erally, the unnormalized state including the pre-factor, 
exp[<&g(t)]|gi, .., qM)\oi q {t)), independently of the pres- 
ence of the quantum jumps provides us a significant 
mathematical simplification. Moreover, we will use the 
result that after the time t > 1/k, all light amplitudes 
a q are constant for all Fock states \q±, .., gjvf)- 

Let us now consider the full initial state given by 
the superposition As stated, the evolution of each 
term is independent and contains the continuous part 
exp[$ g (t)]|gi, ..,q M )\a q {t)). Thus, applying the jump op- 
erators at the times of the photodetections ti,t2,..t m 
leads to the following analytical solution for the condi- 
tional quantum state at the time t after m photocounts: 

|* c (m,t)) = -pTK^2a q (t 1 )a q [t 2 )...a q (t m ) 

\ ' q 

xe^c° q \ qi ,...,q M )\a q (t)), (6) 



where 



F(t) = 




a 9 (ti)| 2 |a 9 (i 2 )| 2 ...|a g (i m )P|e*,(*)p|c0p 



is the normalization coefficient. 

In contrast to a single atomic Fock state, the solution 
([6]) , in general, is not factorizable into a product of the 
atomic and light states. Thus, in general, the light and 
matter are entangled. Moreover, in contrast to a sin- 
gle Fock state, the quantum jump (applying a\) changes 
the state, and the evolution of the full \^> c (m,t)) is not 
continuous. 

The general solution ^ valid for all times simplifies 
significantly for t > when all a q are constants, and 
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if the first photocount occurred at t\ > 1/k, when all a q 
has already become constants. The latter assumption is 
especially probable for the small cavity photon number, 
since the probability of the photon escape within the time 
interval t is 2Kt{a\a\) c . This solution takes the form 

\*c(m,t)) = ^^ a ™ e *' (t)c °l9i.-.9A/>K), (7) 

q 

_ 77 - iU 10 a D q 10 

* g (t) = -\a q \ 2 Kt + {^a* - iU w a a D q w a* q - c.c.)t/2, (9) 
with the normalization coefficient 



F W = ^Z\a q \ 2m e- 2 \ a ^ Kt \c° q \ 2 . 

The solution ([7]) does not depend on the photocount 
times t\,t2, ..t m any more. Note however, that even for 
t > 1/k, when all a q reached their steady states and are 
constants, the solution ([7J is still time-dependent. Thus, 
the time t = 1/k is not a characteristic time scale for the 
steady state of the full solution ([6]) and ([Jj . The station- 
ary light amplitudes a q in (J5J) are given by the Lorentz 
functions in the absolute correspondence with classical 
optics. The function § q (t) has also simplified and con- 
tains the first real term responsible for the amplitudes 
of the coefficients in the quantum superposition, and the 
second imaginary term responsible for their phases. 

The solutions © and ([7]) show, how the probability 
to find the atomic Fock state ..,qm) (corresponding 
to the classical configuration q) changes in time. Such a 
change in the atomic quantum state appears essentially 
due to the measurement of photons and is a direct con- 
sequence of the light-matter entanglement: according to 
quantum mechanics, by measuring one of the entangled 
subsystem (light) one also affects the state of another 
subsystem (atoms). Now we can focus on that purely 
measurement-base dynamics, since other obvious sources 
of the time-evolution (e.g. tunneling) we neglected in 
our model. The initial probability to find the Fock state 
\qx, ..,?m) is Pq(0) = I c° 1 2 . From Eq. ^ the time evolu- 
tion of this probability is given by 

p q (m,t) = \a q (t 1 )\ 2 \a q (t 2 )\ 2 ...\a q (t m )\ 2 (10) 
x\e*«W\ 2 p q (0)/F\t). 

For t,h > 1/k [cf. Eq. 0] it reduces to 

p g (m, t) = \a q \ 2m e- 2 ^ 2 ^ Pq {Q)/F 2 {t). (11) 

In the following, we will demonstrate the applications 
of the general solutions © and J7]) in the examples, 
where, for simplicity, only a single statistical quantity is 
important, instead of the whole set of all possible config- 
urations q. As particular examples, we will consider the 
cases, where that statistical quantity (let us now call it z) 



is the atom number at K lattice sites or the atom num- 
ber difference between odd and even sites. Thus, instead 
of the huge number of detailed probabilities p q (m, t), we 
will be interested in the probability p(z, m, t) to find a 
particular value of z. For the initial state ([2]) at t = 0, 
Po(z) = ^2 q i l c g'| 2 with the summation over all configu- 
rations q' having the same z. As under our assumptions 
all light amplitudes a q (t) depend only on z, we change 
their subscript to z and write the probability to find a 
particular value of z at time t after m photocounts: 

p(z,m,t) = \a z (t 1 )\ 2 \a z (t 2 )\ 2 ...\a z (t rn )\ 2 (12) 
xe 2 ^^p (z)/F 2 (t), 

For t,t\ > 1/k it reduces to 

p{z, m, t) = \a z \ 2m e~ 2 ^ 2Kt Po (z)/F 2 (t), (13) 
F(t) = ^K| 2m e- 2 KI 2Kt p (z). 

In the following we will consider the solution iff]). 
When the time progresses, both m and t increase with 
an essentially probabilistic relation between them. The 
Quantum Monte Carlo method [HJ establishes such a 
relation, thus giving a quantum trajectory. Note, that 
thanks to the simple analytical solution ([7]), the method 
gets extremely simple. The evolution is split into small 
time intervals 5ti. In each time step, the conditional 
photon number is calculated in the state Eq. Q, and 
the probability of the photocount within this time in- 
terval 2K{a\a\) c Sti is compared with a random number 
< e, < 1 generated in advance, thus, deciding whether 
the detection (if 2K{a\a\) c 5ti > ej) or no-count process 
(otherwise) has happened. 



IV. TRANSVERSE PROBING 

In this section we will consider a case, where only the 
transverse probe clq is present, while the probe through 
the mirror does not exist, 7 7 = 0. We presented this situa- 
tion in the previous letter [261 ] . Here we remind the most 
important results and underline the difficulties in the 
state preparation using the transverse probing. In Sec. 
VI, we will switch to a new geometry, probing through a 
mirror, and will show that such a geometry contains the 
features similar to the transverse probing, but is more 
flexible and enables us to solve the difficulties associated 
with the transverse probing scheme. 

In this section we neglect the dispersive frequency shift 
assuming that j7nf n «koi A p . Thus, the light am- 
plitudes a q will only depend on the quantity D\ n . In 
Sec. VI, in contrast, we will focus on the case, where the 
dispersive mode shift is very important. 
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A. Preparation of the atom-number squeezed 
states 

The measurement of photons scattered in the direc- 
tion of a diffraction maximum (Bragg angle) leads to a 
preparation of a state with the reduced (squeezed) fluc- 
tuations of the atom number at the lattice region with 
K illuminated sites [26j]. The condition of the diffrac- 
tion maximum for the scattering of light from the probe 
wave ao into the cavity mode a\ is the following: the 
atoms at all lattice sites scatter the light in phase with 
each other. For the plain standing or traveling waves, 
this condition means that in the expression for the op- 
erator D w = £\ =1 ul(Tj)uo(rj)nj, u*(rj)uo(rj) = 1 for 

all sites j. Thus, Z) 10 = Nk is reduced to the operator of 
the atom number at K illuminated sites. In Eq. ((7|), after 
neglecting the dispersive frequency shift, the only statis- 
tical quantity is giving the atom number at K sites 
for the configuration q. We will call this single statistical 



quantity as z: D\ 



which varies between and N 



reflecting all possible realizations of the atom number at 
K sites. 

From Eq. the light amplitudes in the diffraction 
maximum are proportional to the atom number z: 



a z = Cz, with C = 



iUwao 
(iA p - k) ' 



(14) 



Thus, the probability to find the atom number z is given 
from Eq. ||T3J) by 



p(z,m,t) = z' m e- z T Po (z)/F 2 , 



(15) 



where we introduced the dimensionless time r = 
2|C| 2 k£ and new normalization coefficient F such that 

If the initial atom number z at K sites is uncertain, 
Po(z) is broad. For the SF state the probability to find 
the atom number z at the lattice region of K sites is given 
by the binomial distribution [24| 



N\ ( K\ z ( K 
PSF{Z) = zl(N - z)l [m) 



N-z 



(16) 



For a lattice with the large atom and site numbers 
N, M 3> 1, but finite N/M, it can be approximated as a 
Gaussian distribution 



PSf{z) 



('-'or 



2na z 



(17) 



NK/M and 



with the mean atom number (Nk) = zq 
a z = y/N(K/M)(l - K/M) giving the full width at a 
half maximum (FWHM) 2cr 2 V21n2. The atom number 
variance in the SF state is (AiV^) 2 = (N K ) - (N K ) 2 = 

Eq. (fT5|) shows how the initial distribution po(z) 



its maximum at z\ — y/m/r and the FWHM Sz « 
y/2 In 2/t (for Sz -C Z\), Thus, multiplying po(z) by 
this function will shrink the distribution p(z,m, t) to a 
narrow peak at z% with the width decreasing in time. 

Physically, this describes the projection of the atomic 
quantum state to a final state with the squeezed atom 
number at K sites (a Fock states \z±,N — z\) with the 
precisely known atoms at K sites Z\ and N — z\ atoms 
at M — K sites). Thus, the final quantum state of the 
light-matter system is 



|* c ) = \z!,N -Z!)\a Zl ), 



(18) 



which is a product state showing that the light and mat- 
ter get disentangled. A priori z\ is unpredictable. How- 
ever, measuring the photon number m and time t, one 
can determine z\ of the quantum trajectory. 

Even the final Fock state (fT5|) can contain the signif- 
icant atom-atom entanglement, as this is still a many- 
body state. In general, the Fock state \zi,N — Z\) cannot 
be factorized into the product of two states for K illu- 
minated and M — K unilluminated lattice sites. Thus, 
the entanglement can survive even between two lattice 
regions, which depends on the value of z\ realized at 
a particular quantum trajectory. For some cases, the 
factorization is possible. For example (cf. Ref. [2q). 
the initial superfluid state after the measurement ap- 
proaches the product of two superfluids: \z\, N — Z\) = 
\SF) Z1i k\SF) n _ z1} m_k- 

Note, that our model does not specify how K sites were 
selected, which is determined by the lattice and light ge- 
ometry. The simplest case is to illuminate a continuous 
region. However, one can also illuminate each second site 
by choosing the probe wavelength twice as lattice period 
and get an atom number squeezing at odd and even sites. 
In this way, one gets a measurement-prepared product of 
two SFs "loaded" at sites one by one (e.g. atoms at odd 
sites belong to one SF, while at even sites to another). 
While the initial SF, as usual, shows the long-range co- 
herence (b\bj) with the lattice period, the measurement- 
prepared state will demonstrate the doubled period in 
(b\bj) (bj is the atom annihilation operator such that 

b^bj = hj). Thus, even though our model does not in- 
clude the matter filed operators bj, but only the atom 
number operators hj, the matter coherence can be still 
affected and modified by our QND measurement scheme 
in a nontrivial way. 

The conditioned cavity photon number (a\a±) c (t) — 
\C\ 2 5^2=0 z2 p( z i m i t) i s given by the second moment of 



p(z,m,t). Finally, it reduces to 



\a\ai) 



\C\ 



2 z 2 



rc- 



changes in time. The function 



,2m 



exp (- 



has 



fleeting a direct correspondence between the final atom 
number and cavity photon number, which is useful for 
the experimental measurements. 

Further details of the light measurement at a diffrac- 
tion maximum were presented in Ref. [26| . 
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B. Preparation of the Schrodinger cat states 

The macroscopic superposition state (Schrodinger cat 
state) can be prepared detecting light at the direction of 
a diffraction minimum (2o| . The condition of the diffrac- 
tion minimum for the scattering of light from the probe 
wave do into the cavity mode a\ is the following: the 
atoms at the neighboring lattice sites scatter the light 
with the phase difference tt. For the plain standing or 
traveling waves, this condition means that in the ex- 
pression for the operator Z) 10 = $^7=1 u i( r j) u o( r j)f l ji 



uifojuofo) = Thus, d w = 

is the operator of atom number difference between odd 
and even sites (in this subsection, we consider all sites 
illuminated, K — M). Similarly to the previous subsec- 
tion, in Eq. , after neglecting the dispersive frequency 
shift, the only statistical quantity is Dl Q giving the atom 
number difference for the configuration q. We will call 
this single statistical quantity as z: D\ = z, which varies 
between — N and N with a step 2 reflecting all possible 
realizations of the atom number difference. 

Equations and (fT5|) keep their form for the diffrac- 
tion minimum as well, however, with a different meaning 
of the statistical variable z, which is now a realization of 
the atom number difference, and p{z,m, t) is its proba- 
bility. 

For the SF state the probability to find the atom num- 
ber at odd (or even) sites z [z — (z + N)/2 because the 
atom number difference is z and the total atom number 
is N] is given by the binomial distribution [2~ij 



Psf{z) 



z\(N 



- sv 



Q_ 

M 



N-i 



(19) 



where Q is the number of odd (or even) sites. For even 
M, Q = M/2 and Eq. (0 simplifies. For a lattice with 
the large atom and site numbers N, M 3> 1 , but finite 
N/M, this binomial distribution, similarly to the previ- 
ous subsection, can be approximated by a Gaussian func- 
tion. Changing the variable as z — 2z — N we obtain the 
Gaussian function for the probability to find the atom 
number difference z: 



Psf(z) 



1 



27T(7 Z 



with the zero mean z and a. 



(20) 



N giving the FWHM 



2cr 2 -\/2 In 2. The variance of the atom number difference 
in the SF state is <j 2 z = N. 

The striking difference from the diffraction maximum 
is that our measurement and the probability (fT5|) are not 
sensitive to the sign of z, while the amplitudes a z = 
Cz are. So, the final state obtained from Eq. ([7} is a 
superposition of two Fock states with z\^ = ±-y/m/r 
and different light amplitudes: a Z2 = —a Zl , 

|¥ c ) - ^{W)\a Zl ) + (-in -z x )\- a Zl )). (21) 



In contrast to a maximum, even in the final state, the 
light and matter are not disentangled. In principle, one 
can disentangle light and matter by switching off the 
probe and counting all leaking photons. Then both \a Zl ) 
and | —a Zl ) will go to the vacuum |0). Thus one can pre- 
pare a quantum superposition of two macroscopic atomic 
states (\zi) + (— l) m | — Zi))/V2, which is a Schrodinger 
cat state that, in the notation of odd and even sites, reads 



71 0- 



zi N-zi 



>+(-in 



N — z\ N + zi, 



We have shown that the detection of photons at the di- 
rection of a diffraction minimum leads to the preparation 
of the Schrodinger cat state (f2"Tj) . The physical reason for 
this is that the quantum measurement of photons deter- 
mines the absolute value of the atom number difference 
\z\. However, as the photon number is not sensitive to 
the sign of z, one ends up in the superposition of states 
with the positive and negative values of z. 

Unfortunately, Eq. (|2ip demonstrates a very strong 
disadvantage of such a method, which makes it difficult 
to realize experimentally. Each photodetection flips the 
sign between two components of the quantum superpo- 
sition in Eq. (f2Tj) . This means, that if one loses even a 
single photocount, which is very probable for a realistic 
photodetector, one ends up in the mixture of two state 
(I2ip with plus and minus signs. Such a mixed state does 
not contain any atomic entanglement any more, in con- 
trast to the pure state (|2"Tj) , which is a highly entangled 
one. 

Formally, the appearance of the sign flip (— l) m in 
Eq. (l2"Tj) originates from Eq. (O, which contains the co- 
efficient a™ in the pre-factor of each Fock state. Since, 
in the final state, two components with opposite signs 
of light amplitudes survive (a Z2 = —a Zl ), the term a™ 
produces the coefficient (— l) m . Therefore, even one pho- 
todetection changes the phase between two components 
in Eq. (|2ip by Aip± — tt, which is the maximal possi- 
ble phase difference. The idea to make the preparation 
scheme more stable with respect to the photon losses is 
based on the possibility to make this phase jump A<pi 
smaller than tt. In this case, loosing one photon will also 
lead to the mixture of two cat states. However, if those 
cat states are not very different (i.e. the phase difference 
Aipi is small), the mixed state will still contain significant 
atomic entanglement. 

In Sec. VI, we present a scheme to prepare the 
Schrodinger cat state with a phase difference between 
two components smaller than tt. Thus, such a scheme 
is more practical and robust with respect to the photon 
losses and, hence, decoherence. 

The atom number squeezed states (|18p prepared by 
observing light at a diffraction maximum is indeed more 
robust than the Schrodinger cat state (|2"Tj) obtained at 
a diffraction minimum, as the former do not have any 
phase jump. However, the convenient property of the 
measurement at a minimum is that during the same 
time interval (e.g., during the shrinking time which is the 
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same for a maximum and minimum, Sz ~ y^2 In 2/r) the 
number of photons scattered at a diffraction minimum 
((a\ai) = \C\ 2 N) is much smaller than the one scattered 
at a maximum ((ojoi) = \C\ 2 N%) (H [H, SH] . 

Further details of the light measurement at a diffrac- 
tion minimum were presented in Ref. (26[. 

V. PHOTON STATISTICS 

In this section, we consider three kinds of photon 
statistics: (i) statistics p$(n,m,t) of the photon number 
n in a cavity after m photons were detected outside the 
cavity, (ii) statistics P{m, t) of the photocount number 
m, and (iii) statistics Pr(m,t) of the photocount num- 
ber m if, after the time measurement T, photons 
were detected. 

First, let us consider the statistics of the number of 
photons in a cavity after the measurement time t and m 
photodetections. The joint probability to find a number 
of photons in a cavity n together with finding the atomic 
state in the Fock state \q) = \qi, .., qm) is obtained by 
projecting the general solution ([6|) on the state \q}\n} 
and is given by 

W q (n,m,t) = M^e-I^WI^Kt), 

where the probability to find the atomic Fock state 
p q (m,t) is given by Eq. (fit)]) . However, the probabil- 
ity to find n photons in a cavity independently of the 
atomic state is obtained by projecting the solution (J5|) 
on the light Fock state \n) and taking the trace over the 
atomic states. Thus, the cavity photon number distribu- 
tion function is given by 

»(u, ffl ,()^wf: e -w<)i ! p , (mii) , 

z — ' n\ 

g 

where the sum is taken over all possible configurations q. 

If, as in the previous section, the only atomic statistical 
quantity is z, the sum is simplified and the probability 
p$(n,m,t) to find n photons in a cavity after the mea- 
surement time t and m photodetections reads as 

Wn, m,t)=y M*)l 2 V l"*WI 2 p( Zi m, t), (22) 

z 

where the atom number distribution p(z, m, t) is given by 
Eq. CEl- 

In general, p$(n, m, t) is a super-Poissonian distribu- 
tion. During the measurement, the atomic distribution 
p(z,m,t) shrinks to one or two symmetric peaks corre- 
sponding to the atom-number squeezed or macroscopic 
superposition states. Thus, after some time, only a single 
term (or two equal terms) survives in the sum, and the 
cavity photon statistics p$(n, m,t) evolves from super- 
Poissonian to Poissonian one. This fact can be checked 
experimentally. 



During the measurement, the mean conditional pho- 
ton number (a\ai) c approaches the value |a Zl | 2 , which 
enables one to determine the final atom number z\ by 
measuring the photon number in a cavity. 

Let us now consider the probability P(m, t) to detect 
m photons within the time t, if the initial atom number 
distribution is Po{z). As shown, for example, in Ref. [5^ |. 
this probability can be obtained from the state ([6]) us- 
ing the integration over all detection moments tx,t2,-,t m 
from to t (because we are not interested in the time mo- 
ments, but only in the total number of the photocounts 
m) and taking the trace. The simple result can be ob- 
tained for the case t, t\ > 1/k, since the solution (J7J) docs 
not depend on the detection times. The probability reads 

PKt) = EM^e-l-^o ( z), (23) 
* — ' ml 

z 

where the powers of m appear due to the m time- 
integrations. 

In contrast to the probability p$(n,m,t), Eq. (|2"2"j) . 
which characterizes the conditional distribution of the 
cavity photons at a particular trajectory, the probability 
P(m, t) depends on the initial atom number distribution 
Po(z) and is not a characteristic of a particular quantum 
trajectory, but rather of an ensemble average. In gen- 
eral, this is a super-Poissonian distribution, which does 
not approach any Poissonian one. From Eq. (|23[) . the in- 
crease of the mean photocount number with time is given 
by 

(m) = 2Kt(a\ai) , 

where (a\ai)o is not a conditional photon number, but 
the one calculated for the initial atomic state. In aver- 
age, the photocount number linearly increases in time. 
However, as the distribution is not Poissonian, the fluc- 
tuations of the photocount rate m/t do not decrease to 
zero. 

One can also introduce another statistical distribution: 
the distribution Pr(m, t) of the photocount number m 
if, after the measurement time T, wit photons have been 
already detected. Similar approach as in Eq. (j23|) . leads 
to the following result: 

PT( m ,t)=E [2K|a ' |2(t ," Tr ^^' |i(t - r) 

^—^ ml 

z 

xp(z,m T ,T), (24) 

which, in contrast to Eq. (|23[) . depends not on the ini- 
tial atomic distribution, but on that at time T. Thus, 
this probability combines the quantum trajectory evolu- 
tion up to the time T, and the ensemble average after 
that. As we know, the atomic distribution p(z,m,T,T) 
approaches the single peak for the Fock state with in- 
creasing T. Therefore, this photocount probability ap- 
proaches the Poissonian distribution with increasing T, 
the fluctuations of m grows in time as yfm ~ y/t, and the 
fluctuations of the photocount rate m/t vanishes with in- 
creasing time. 



9 



VI. STATE PREPARATION BY CAVITY 
TRANSMISSION MEASUREMENT 

We now switch to a different probing scheme, using 
the probe through the mirror with the amplitude 77 (Fig. 
1), and assuming no transverse probe, clq = 0. From 
Eq. ([8]) we see that the light amplitudes depend only on 

the single statistical quantity D\ x — J2j=i \ u i( r j)\ 2( lji 
which reduces to the atom number at K sites for the 
traveling wave a\ at any angle to the lattice, or for the 
standing wave a\ with atoms trapped at the antinodes. 
Thus, in this case, the statistical quantity is z = D\ 11 
which changes between and N. The term UnD^ in 
Eq. ([5]) has a meaning of the dispersive frequency shift of 
the cavity mode, due to the presence of atoms in a cavity. 

The light amplitude a z from Eq. ([5]) can be rewritten 
as a function of the atom number at K sites z: 



c 



— — ttt~ with C' = — , 

1{Z — Z p ) + K/Ull K 



(25) 



and the parameter z p = A p /Un is fixed by the probe- 
cavity detuning A p . The probability to find the atom 
number z (11311 takes the form 



e -r' (k/Uh) 2 /[{z-z p ) 2 + {k/Uh) 2 } 

P(^m,t) = [{z _ Zp) 2 + {K/Ull) 2 ]m P^)/F' 2 l26) 

where the dimensionless time is r' ~ 2\C'\ 2 Kt. 

The pre-factor function in front of po(z) in Eq. (j2l)|) 
has a form more complicated than the one we had for 
the transverse probing in Eq. (|15[) . However, it provides 
us a richer physical picture. In contrast to the transverse 
probing, this function allows us a collapse to both the 
singlet and doublet distribution. Thus, the measurement 
of the cavity transmission contains both cases of trans- 
verse probing at a diffraction maximum and minimum 
considered in Sec. IV. 

If for a particular quantum trajectory the number of 
photocounts is large such that m/r 1 > 1, the distribution 
p(z,m,t) collapses to a single-peak function at z p . This 
singlet shrinks (FWHM) as 5z « 2{k/Uu) f/2 In 2/r', 
which can be very fast for a high-Q cavity with the small 
k. In the estimation of 8z, we used the assumption that 
the peak has already become narrow, Sz <C k/Uh, and 
that to « 2K(a[ai) c t. Similar to the transverse probing, 
in the final state, the light and atoms are disentangled: 
\z p , N — z p )\a Zp ) , which corresponds to the atom-number 
squeezed (Fock) state and coherent light state. 

If, however, the number of photocounts is small, 
m/r' < 1, the distribution p(z,m,t) collapses to a dou- 
blet centered at z p with two satellites at zi_2 — z p ± Az 
with 



K IT „ 

Az = —J--1. 
U11 V m 



(27) 



Az), each of its component shrinks in time 



as 



Az 1 



/21n2 



C/^Az 2 , 



U^Az- 



Physically, tuning the probe at A p , we may expect 
scattering from the atom number z p providing such a 
frequency shift A p . If the photocount number is large 
{m/r' > 1), indeed, the atom number is around z p and 
it collapses to this value. However, if m is small, we 
gain knowledge that the atom number z is inconsistent 
with this choice of A p , but two possibilities z < z p or 
z > z p are indistinguishable. This collapses the state to a 
superposition of two Fock states with Z\£, symmetrically 
placed around z p . 

Thus, the transmission measurement scheme allows 
one to prepare both the atom number squeezed state and 
the Schrodinger cat state. The appearance of the singlet 
for squeezed state or the doublet for cat state can be 
determined by measuring m and i, or the final photon 
number (a\ai) c = \a Zl \ 2 . From Eq. ([23)) . we see that 
there is a direct correspondence between the final cavity 
photon number and the parameter of the doublet Az: 



{a\ai) c = 



U^Az 2 



(28) 



When the doublet has become well-separated (i.e. 5z <C 



which makes the experimental determination of the dou- 
blet position possible. The parameter Az can be also 
determined by measuring the photocount number m and 
time t using Eq. (j2"7) . 

In Figs. 2-5, we present the results, where the quan- 
tum trajectories of qualitatively different kinds were re- 
alized. In all figures, the initial state is the superfluid 
with the atom number N = 100 at M = 100 lattice 
sites, the half of all lattice sites K = 50 are illumi- 
nated by the cavity mode. The initial distribution of 
the atom number at K sites is given by Eq. (|16p and can 
be well approximated by the Gaussian distribution Q17p 
with the mean value (Nk) = zq = NK/M = 50 and 
<r z = y/N(K/M)(l- K/M) = 5. 

Figure 2 presents the results for a quantum trajec- 
tory, which leads to the collapse to a single-peak dis- 
tribution. In this case, the probe is detuned in such way, 
that the probe-cavity detuning corresponds to the center 
of the atom number distribution: z p — A p /U\i — 50. In 
Fig. 2(a), the evolution of the atom number probability 
p(z,m,t), Eq. (|2l))) . is shown. The curve A is the initial 
atom number distribution having the Gaussian shape. 
Curve B shows the non-Hermitian evolution of the prob- 
ability just before the first jump (m — 0). As, at this time 
interval, the number of photocounts is small (m/r' < 1), 
the distribution tends to a doublet. This is assured by 
the exponential factor in Eq. ([23)) leading to the suppres- 
sion of the component at z = z p as one does not record 
the photocounts at the expected detuning. However, just 
after the photocount (to = 1) this distribution instantly 
changes to the curve C, which has already a single peak 
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FIG. 2: (Color online) Photodetection trajectory leading to a single-peak distribution (atom number squeezing). The probe- 
cavity detuning is chosen such that z v — 50 coincides with the initial distribution center, (a) Shrinking atom number distribution 
at different times r' = (number of photocounts is m = 0), 0.7 (just before the first count, m = 0), 0.7 (just after the first 
count, m — 1), 1.1 (m = 1), 1.1 (m = 2), 14.6 (m = 17) (A-F); (b) time dependence of the photocount number m, the red 
fine is for m = r'; (c) decreasing width of the atom number distribution; (d) reduced conditioned photon number (a\a\) c /\C'\ 2 
with quantum jumps; (e) photon Mandel parameter. Initial state: SF, N = 100 atoms, M = 100 lattice sites, K — M/2 — 50 
illuminated sites. 



as for this trajectory it turns out that m/r' > 1. The 
switch to a single peak is assured by the Lorentzian fac- 
tor in Eq. ([26)) . After that, before the second jump, the 
probability decreases at z = z p and broadens (curve D), 
but jumps upwards again and narrows, when the second 
jump occurs (curve E). Finally, after many jumps, the 
probability distribution becomes narrow and has a single 
peak at z = z p (curve F). 

Figure 2(b) shows the number of photocounts m grow- 
ing in time. It is clear that m stays always near the line 
m = t' . The appearance of the singlet is assured by the 
fact, that at the initial stage m/r' > 1. 

Figure 2(c) shows the evolution of the width of the 
atom number distribution, yj (AN K ) 2 , which decreases 
to zero reflecting the shrinking distribution. Figure 2(d) 
shows the reduced conditioned photon number in the cav- 



ity (a\ai) c /\C'\ 2 . One sees that for the initial atom dis- 
tribution it starts from a relatively small value. However, 
as the atomic state goes to the Fock state with atom num- 
ber z p exactly corresponding to the detuning A p , it ap- 
proaches the maximal possible value (a\ai) c /\C'\ 2 = 1. 
Moreover, one can easily see the quantum jumps in 
the initial stage, when the light field consists of several 
coherent-state components. Finally, the jumps disappear 
as the light state approaches a single coherent state \a Zp ), 
when the atomic state approaches a Fock state. It is in- 
teresting to note, that the photon escape from the cav- 
ity (photocount) leads to the increase of the conditional 
photon number in the cavity, while the no-count process 
leads to its decrease. This is a counter-intuitive char- 
acteristic feature of the super-Poissonian photon statis- 
tics and is determined by the conditional nature of the 
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FIG. 3: (Color online) Photodetection trajectory leading to a doublet distribution (Schrodinger cat state). The probe-cavity 
detuning is chosen such that z v — 50 coincides with the initial distribution center, (a) Shrinking atom number distribution at 
different times r' = (number of photocounts is m = 0), 16.4 (just before the first count, m = 0), 16.4 (just after the first 
count, m = 1), 58.2 (m = 1), 58.2 (m = 2), 2017.6 (m = 73) (A-F); (b) time dependence of the photocount number m, the red 
line is for m — r'; (c) decreasing width of one of the peaks in the atom number distribution; (d) reduced conditioned photon 
number (a\a\) c /\C'\ 2 with quantum jumps; (e) photon Mandel parameter. Initial state: SF, N = 100 atoms, M = 100 lattice 
sites, K — M/2 = 50 illuminated sites. 



probabilities considered [59j. Figure 2(e) shows the re- 
duced Mandel parameter Q/\C'\ 2 characterizing the pho- 
ton number variance: Q — ((n|,) — (n$) 2 )/(n$) — 1, where 
n<p = (a\ai) c is the conditioned photon number. One can 
see how it decreases to zero corresponding to the coherent 
state of light. The parameters shown in Figs. 2(b)-2(e) 
can be measured experimentally thus presenting the ver- 
ification of our theory. 

Figure 3 presents the results for the case, where the 
detuning also corresponds to the atom number distribu- 
tion center z p = 50, but the quantum trajectory leads to 
the doublet distribution (Schodinger cat state). Figure 
3(a) shows the evolution of the atom number distribution 
p{z, m, t). Curve A is the initial Gaussian distribution. 
Before the first jump (curve B) the distribution evolves 
to a doublet-like, similar to the curve B in Fig. 2(a). 



However, in contrast to Fig. 2(a), the first photocount 
does not return the distribution back to the single peak, 
and it stays doublet-like (curve C). This is so, because 
the first jump occurs rather late such that m/r' < 1. 
Next jumps are rather late as well, so the distribution 
evolves to a doublet (curves D, E, and F). 

Figure 3(b) shows the number of photocounts m grow- 
ing in time. It is clear that, in contrast to Fig. 2(b), m 
is always much smaller than the line m = r' (red line), 
which gives the experimental possibility to claim the ap- 
pearance of the doublet. 

Figures 3(c), 3(d), and 3(e), similarly to Fig. 2, show 
the decrease of the width of one of two peaks of the atom 
number distribution, conditioned cavity photon number 
with disappearing jumps, and the Mandel parameter ap- 
proaching zero. Note, that in contrast to Fig. 2(d), the 
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FIG. 4: (Color online) Photodetection trajectory leading to 
a single-peak distribution (atom number squeezing). The 
probe-cavity detuning is chosen such that z v — 60 is at the 
wing of the initial distribution. Shrinking atom number distri- 
bution at different times, (a) r' = (number of photocounts 
is m — 0), 0.7 (just before the first count, m = 0), 0.7 (just 
after the first count, m = 1) (A-C); (b) t' = (m = 0), 1.1 
(m = 1), 1.1 (m = 2), 14.6 (m = 17) (A-D). Initial state: 
SF, N = 100 atoms, M = 100 lattice sites, K = M/2 = 50 
illuminated sites. 



conditioned photon number (a[ai) c /\C'\ 2 does not ap- 
proach the maximal value 1, but rather decreases to a 
smaller value given by the doublet splitting Az, Eq. (|2"5)l . 
Thus, the appearance of the doublet can be characterized 
by experimentally measuring m and t [Fig. 3(b)] or the 
cavity photon number [Fig. 3(d)]. 

Figures 4 and 5 present another situation, where the 
probe-cavity detuning is chosen such that z p = 60, which 
corresponds to a wing of the atomic distribution function. 

Figure 4 shows the evolution of the atom number distri- 
bution in case, where it collapses to a singlet at z p . Thus, 
here even the conditioned mean atom number changes 
from 50 to z p = 60. Other characteristics of this process 
look very similar to the ones presented in Figs. 2(b)-(e). 



Figure 5 shows the collapse to a doublet distribution 
around z p — 60 with Az — 7. Thus, two satellites in 
Fig. 5(a) are placed at z\ — 67 and Z2 = 53. However, 
while the satellite at z 2 = 53 is near the maximum of the 
initial distribution and is well seen, the second satellite at 
zi = 67 falls on the far wing and is practically invisible. 
As a result, the final distribution looks as a singlet at 
Z2 = 53, while the second satellite is very small. The 
fact that one has indeed a doublet can be verified by 
measuring the photocount number [Fig. 5(b)], which is 
obviously less than r', or by measuring the cavity photon 
number [Fig. 5(e)], which is less than the maximal value 
1 and depends on Az. The measurement of the mean 
atom number [Fig. 5(c)] or width of atomic distribution 
[Fig. 5(d)] would not distinguish between the singlet and 
doublet. 

This suggests us a method to prepare the macroscopic 
superposition state with unequal amplitudes of two com- 
ponents. Choosing the detuning A p such that z p is not 
at the center of the initial atom number distribution, one 
can expect that one of the satellites will be very probable 
near the distribution center, while another one will fall 
on its wing. 

The most important difference of the Schrodinger cat 
state prepared by the transmission measurement from 
the one prepared by transverse probing (|21[) is that the 
phase difference between two components of the quantum 
superposition is not limited by the values of and 7r as in 
Eq. (|2ip . Although the light amplitudes corresponding 
to z\ and zi has equal absolute values, their phases are 
opposite: a Zl = \a Zl | exp (iip) and a Z2 = \a Zl | exp (— iip), 
where from Eq. (j2"5|) the light phase cp is 



ip = — arctan 



UnAz 



(29) 



Moreover, the imaginary parts of the functions $ q in 
Eqs. ([7]) and ^ are also different for z\ and Z2 and 
have opposite signs: <&(<) = Im$ Zl (i) = — Im$ 22 (i) = 
Im(rya*J. Using Eq. ([2>]), 

$(<) = {a^UuAzt. (30) 
Thus, using Eq. 0, the cat state is given by 

\^c) = ^[e m ^ + ^\z 1 )\a Zl )VM^) 



+e- imv - im \z2)\a Z2 )VM^)}, (31) 

which is a macroscopic superposition of two Fock states 
with the atom numbers zi and z 2 at K sites. 

Using the probe-cavity detuning A p one can chose the 



center of the doublet z p = A p /Uu. Moreover, Eq. |3T 
shows that using this detuning one can influence, at least, 
probabilistically, the ratio between two components in 
the superposition (f3Tj) . In particular, if we tune the 
probe frequency such that z p coincides with the center 
of the initial atom number distribution po(z), the prob- 
abilities of the symmetric doublet components will be 
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Time t' Time r' 

FIG. 5: (Color online) Photodetection trajectory leading to a doublet distribution (Schrodinger cat state). The probe-cavity 
detuning is chosen such that z v = 60 is at the wing of the initial distribution, (a) Shrinking atom number distribution at 
different times r' = (number of photocounts is m = 0), 5.4 (m = 1), 163.5 (m = 3), 2006.9 (m = 39) (A-D); the arrow shows 
the position of the very small doublet component, which is almost invisible; (b) time dependence of the photocount number 
m, the red line is for m = r'; (c) conditioned mean atom number; (d) decreasing width of the atom number distribution; (e) 
reduced conditioned photon number (a\a\) c /\C'\ 2 with quantum jumps. Initial state: SF, N = 100 atoms, M = 100 lattice 
sites, K — M/2 = 50 illuminated sites. 



equal, po(zi) = pq(z2), providing the equal amplitudes of 
the cat components in Eq. ([3"T]) . 

The phases $(i) and mip have opposite evolution in 
time. The pase term Eq. (|30[) . grows in time linearly 
and dctcrministically, while the term mip grows in time 
stochastically according to the growth of the photodetec- 
tion number to. In average, (to) = 2n\a Zl \ 2 t. Thus, in 
general, the phase difference between two cat components 
grows in time linearly. However, for some parameters, 
the growth of the average (m)ip and &(t) can be com- 
pensated. From Eqs. and (m)<p + $(t) = for 
UhAz/k » 2.33. So, for the particular cat components, 
the phase difference does not grow in average. However, 
as to is a stochastic quantity, its uncertainty grows in 
time as y/t (cf. Sec. V). Thus, the phase difference will 
still grow in time as y/t, which is however much slower 



than the linear growth. 



The problem of photon losses addressed in the previ- 
ous section is related to the stochastic quantity mtp. A 
single photocount changes the phase difference between 
two components by Aipi — 2(p. In contrast to the state 
(|21[) . where the phase jump is always maximal Aip± = n, 
here, this jump can be rather small, if the condition 
UhAz/k < 1 is fulfilled. This means, that to provide 
the robustness of the state with respect to the photon 
losses, the doublet should not be split too strongly. In 
the next section we quantitatively analyze the robustness 
of the cat state. 



14 



VII. ROBUSTNESS OF THE SCHRODINGER 
CAT STATES 

The macroscopic superposition state ([51]) is a pure 
state. In principle, if the measurement is perfect and 
all photons m leaking the cavity are counted by a pho- 
todetector, it will evolve according to Eq. (|3"Tj) staying 
pure. However, if one loses one or more photons, the 
state becomes a mixture of several states corresponding 
to several lost counts I. Thus, if L photons are lost, the 
state is a mixture of L + 1 states of the following form 
for < I < L: 



I*, 



^+^\z 1 )\a Zl )+e~ a ^\z 2 )\a Z2 )}, (32) 



where, for simplicity, we assumed the symmetric super- 
position with po(zi) — Po{z2) and included all known 
phases [for m measured photons and deterministic 
in to the term 7 = mip + <fr(t). The density matrix of the 
state 




1 2 3 

Doublet splitting U^Az/ k 

FIG. 6: (Color online) Purity of the mixed state as a function 
of the doublet splitting Az for different numbers L of the 
photocounts lost: L — 0, 1,3, and 10. 



Pi = ^(\ z i)\^z 1 )(zi\(a Zl \ + \z 2 )\a Z2 )(z 2 \(a Z2 \ 

+e m ^\ Zl )\a Zl )(z 2 \(a Z2 \ 
+e-^-^\z 2 )\a Z2 )(z 1 \(a Zl \). (33) 

The density matrix of the mixture state describing L lost 
photons is given by a sum of the density matrices (I33[) : 

P (L) = ~{\zi)Wzx){zi\{uzA + \z 2 )\a Z2 )(z 2 \(a Z2 \ 



J2j 



L 



1=0 



e mv \ Zl )\a zl ){z 2 \{a Z2 \ 



-lit 



L 



Y Y,^ 2lV \^)\a Z2 )( Zl \(a Zl \). (34) 



1=0 



The quantity characterizing how close is a mixture 
state to a pure state is the so-called purity: P = Tr(p 2 ). 
For a pure state it is maximal and equal to 1, while for 
a maximally mixed state it is minimal and equal to 1/2 
(in our case of the two-component states). The purity of 
the state ([M)) is given by 



Pl = 



1 



1 



1=0 



(35) 



where the sum can be calculated leading to the following 
result for the purity of the mixed state corresponding to 
L lost photons: 



1 

Pi =2 



1 sm 2 (L + l)ip 



(L + l) 



■ 2 

sin tp 



(36) 



For example, in the simplest case, where one photon 
is lost, the density matrix of the mixed state is given by 



the sum of two terms (l33l) with I = and 1 



P (1) = ^(I^KiX^iK^il + \ z 2l\a Z2 ){z 2 \{a Z2 \ 

+ le-^(l + e-^)\z 2 )\a Z2 )(z 1 \(a Zl \), (37) 



which has a purity 
Pi = ^(1 + cos 2 if) 



l + {U n Az/ K ) 2 



(38) 



Equations ([3T)l and (f3"B"]) show that if the phase jump 
associated with the one-photon lost A(pi — 2ip is max- 
imal, Aipi — 7T, the state (|37j) is maximally mixed, be- 
cause all non-diagonal terms responsible for the quantum 
coherence between the states |zi)|a 2l ) and |z2)|a^ 2 ) are 
zero. Its purity (|38p is 1/2, thus no entanglement sur- 
vived after the single photon lost. This is a situation 
of the transverse probing in the diffraction minimum, 
Eq. (|21|) . where the phase jump is 7r, which makes the 
preparation scheme practically difficult. 

However, if ip is small, the purity (|38p and the non- 
diagonal coefficients in Eq. (|3"7| can be rather large and 
close to 1. Thus, after the photon has been lost, one 
gets a mixed state, but that of the high purity. More 
generally, for L photon losses, Eqs. (fM)) and (|3"6")l . the 
total phase jump A(p^ — 2Ltp should be small. Using the 
expression for ip (|29p . one can estimate the condition for 
the high purity of the mixed state as tp < tt/(2L), which 
for a small tp is approximately the condition UnAz/n < 
7r/(2L). As a result, for the doublet, which is split not 
too strongly, the high purity can be preserved even with 
photon losses. 

The purity (|36p as a function of the doublet splitting 
Az for L — 0, 1, 3, and 10 photon losses is shown in 
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Fig. 6, where <p is given by Eq. ((29)) . One can see the 
decrease of the purity with increasing doublet splitting 
and number of photons lost. Note the non-monotonous 
character of its decrease, which means that the larger 
splitting does not automatically leads to a smaller purity. 

VIII. CONCLUSIONS 

We explicitly calculated the quantum state reduction 
induced by measurement of off-resonant scattered light 
from an ultracold quantum gas t rapped in an optical lat- 
tice. In our previous papers [23L[24l25|. the expectation 
values of various atomic and light quantities were ana- 
lyzed, thus, assuming the repeated measurement to en- 
able the averaging procedure. In contrast, in this paper, 
we were focused on a single run of the optical measure- 
ment, i.e., the evolution of the quantities at a single quan- 
tum trajectory. As the scattered light is entangled to the 
atoms, quantum back-action of the light measurement al- 
ters the atomic state. The geometry, which determines 
the form of entanglement, thus dictates the possible mea- 
surement results and final atomic states. 

Any quantum state related to an eigenstate of the 
weighted atom-number operators Di m can thus be pre- 
pared in a probabilistic way. The type of the states can 
be chosen by the optical geometry and their probabilities 
are determined from the initial distribution. Typically 
light detection at the diffraction maximum leads to the 
preparation of the atom-number squeezed states, while 
the detection at a minimum prepares the macroscopic 
superposition states peaked at a pair of atom numbers. 

The robustness of the resulting Schrodinger cat states 
with respect to undetected photons was analyzed. The 
transmission measurement scheme was shown to prepare 
the states more robust than the ones prepared by the 
detection at a diffraction minimum. 

In contrast to recent results in spin squeezing and 
preparation of the spin superposition states, which can be 
also obtained for thermal atoms [13,111,111,13113,1111, in 



our work, quantum nature of ultracold atoms is crucial, 
as we deal with the atom number fluctuations appearing 
due to the derealization of ultracold atoms in space. 

We demonstrated the time evolution of various measur- 
able quantities appearing exclusively due to the measure- 
ment procedure (as other obvious sources of dynamics 
such as tunneling were neglected) . The quantum dynam- 
ics is governed by the quantum jumps and conditional 
evolution. The quantum state preparation is probabilis- 
tic. However, it can be generalized by including the feed- 
back loop, which will enable ones the quasi-deterministic 
state preparation. In this case, the trapping potential 
should be continuously modified depending on the out- 
come of the photodetector. For example, the detection 
of photons at a diffraction maximum squeezes the atomic 
number at K sites around some value z\ , which was not 
known a priori. The potential can be continuously tilted 
in a way to provide the increase or decrease of this atom 
number to enable ones to obtained the number squeezed 
state with a mean value i,\ given a priori. The same 
method can be applied for the atom number squeezing 
at odd or even sites. 

Our results can be applied for other quantum arrays as 
well, e.g., ion strings [55j, and be useful for the prepara- 
tion of various atomic and photonic multipartite entan- 
gled states [60(. Moreover, they can have a relation to 
the problem of interaction of ultracold atoms with other 
bosonic particles, besides photons [6lll62l]. 

Cavity QED with quantum gases can operate with the 
atom numbers ranging from millions to one [63| . Thanks 
to the recent experimental breakthroughs 0, H, [1] , 
preparing various kinds of atom number squeezing is al- 
ready doable, and creation of, at least, Schrodinger "kit- 
tens" [65[ may become practical. 
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